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The following abbreviations are used: 



! - r. 



The coefficients a, 0 , y are evaluated by means of the series (22, . 
23, 24) and the function /(c, k 2 ) is defined in (31). 
One finds 
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Of = -yF c (l) - ydc, + *z)*>c(j) - 4c/f c 
In the above the definition is used 



(36,11) 
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The use of a programmable calculator enables one to compute 
the coefficients in a minute, provided that | k , | and | k 2 | lie below 
0.95. For larger values, the time of computation increases 
rapidly. 



The Limiting Case of equal Cylinder radii R l = R 2 

The formulae simplify considerably in this case. We quote 
them here, but skip the explicit calculations. Although these 
may be derived from the general formulae given above as a 
special case, it is far simpler to start with an integration over 
; m 0, as described in part I. 
For /J, = R 2 , one obtains 



(37,10 



By repeated use of the functional equation given under the 
definition (21), one gets the general relation 



- + (a- D-Ztftfea- 1) 



1 — c 

+ (a - 2) - Zt-^-ttofl - 2) + . . . 

+ (a-*-l)-Z 0>0 (c.a-*-l) 

(A: = 0,1,...) (38.11) 

(33) enables us to transform the series D c (a) and H e (a) into 
those that converge better in the neighbourhood of 1. With a 
programmable calculator, the coefficients or, 0, y can thus be 
determined with great accuracy for values of | c \ up to 0.99 in a 
matter of minutes. 
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Solubility and Diffusion of Chlorine in Silica Glass 

Wilhelm Hermann, Hans Rau, and Joachim Ungelenk 

Philips GmbH Forschungsiaboratorium Aachen, D-5100 Aachen, Germany 

Diffusion / Glasses 

Chlorine diffusion profiles in suica glass have been studied by SEM techniques. From the results the diffusion coefficient D Oog io (0cnr 2 s) = 
-11080 K/T -3.498), the equilibrium solubilitv of CI in silica glass (025 or 0.32 weight per cent at 1 bar Cl 2 pressure and 1438 K or 
1970 K), and the pressure dependence of the sotubiUty were derived. The latter is proportional to the fourth root of the Cl 2 pressure. 
There are no detectable differences in the behaviour of chlorine in silica made by the PCVD method (plasma activated chemical vapour 
deposition) or in commercially produced silica (HERALUX, wave guide quality), but the chlorine concentration is much higher in PCVD 
silica than expected from the small chlorine pressures present during deposition. 



Introduction by Greeuw and Hasper [2] who implanted chlorine in Si0 2 

Very little is known of the solubility and the diffusion of layers grown on silicon. They found that chlorine becomes 

chlorine in silica. In the "Handbook of Glass Data" [1] mobile in Si0 2 above 1000 K and derived a tentative dif- 

nothing is said on this subject. Some information is given fusion coefficient. Hanafusa and Tajima [3] measured 
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chlorine concentration profiles in optical fibres based on 
Si0 2 . 

During deposition of (doped) silica layers in a silica tube 
for the preparation of optical fibres using the PCVD method 
[4] significant amounts of chlorine are incorporated 
(0.3 — 1.1 weight per cent). Since this effect is of rather great 
significance for the PCVD process we decided to determine 
the solubility of chlorine in silica as a function of Cl 2 pres- 
sure and temperature and the diffusion coefficient as a func- 
tion of temperature. 

Experimental 

Samples Prepared by the PCVD Method 

A mixture of oxygen (22- 10"* mol/s) and silicon tetrachloride 
(3.7- 10" 5 mol/s) flowing through a silica tube of 14 mm internal 
diameter at 10 mbar pressure and 1 530 K was activated by a moving 
non-isothermal plasma generated using a microwave resonator 
(2.45 GHz, 500 W). Thus, the reaction 

SiCU + 0 2 - Si0 2 + 2C1 2 (1) 

took place in the plasma region and the Si0 2 formed was con- 
densed at the wall of the tube to clear thin glass layers. After having 
deposited about 03 mm in thickness the reaction was stopped and 
the tube cooled to room temperature. Slices of the tube taken per- 
pendicular to the axis were annealed at different temperatures for 
a sufficiently long time to allow the chlorine incorporated in the 
deposited layers to diffuse into the original silica tube which was 
known to be free of chlorine at the beginning (HERALUX wave 
guide quality). The slices were embedded in suitable plastic and 
polished to optical quality. After vapour deposition of carbon the 
chlorine concentration profiles of the samples were measured with 
a scanning electron microscope (SEM) (runs 1-3). 

Annealing of Silica Ampoules or Tubes with Chlorine Gas 

Known amounts of chlorine were condensed into silica ampoules 
of known volume which were subsequently sealed. Each ampoule 
also contained a rod of silica glass. The ampoules were then an- 
nealed in an electrical furnace. After cooling down and condensing 
the chlorine gas the ampoules were opened and slices were taken 
from the ampoule walls, embedded and polished as described above 
(runs 4 - 6). The rods were treated in a nuclear reactor with 3.2 • 10 1 1 
neutrons cm* 2 s~ l for one hour and the total chlorine diffused 
into them was calculated from the intensity of its characteristic 
gamma radiation (runs 8-10). 

At 1970 K the treatment with chlorine could only be performed 
by using a flow of chlorine gas of exactly ambient pressure through 



a tube mounted in a suitable furnace because silica is no longer 
hard at this temperature (run 7). 

Scanning Electron Microscope (SEM) Analyses 

The experimental conditions in the energy dispersive X-ray anal- 
ysis carried out in a scanning electron microscope were as follows: 
1 2 keV excitation energy, 200 nm spot size, approx. 2 - 1 0 " 9 A probe 
current and 15 degree specimen tilt angle. The method was cali- 
brated by analysing the same specimens with wet chemical methods. 
The profiles were obtained either by point to point measurements 
or, in order to avoid the influence of small defects in the surface 
on the results, by analysing a set of parallel line scans perpendicular 
to the slope of the profile. The propagated 1-sigma statistical error 
achieved was 9% at chlorine levels of 0J weight per cent in silica 
glass. In analysing a chlorine step profile (composite specimen con- 
sisting of a piece of chlorine doped silica glass and a piece of pure 
silica glass precisely polished and fitted to each other) it has been 
verified that under the chosen experimental conditions the spatial 
resolution is better than one micrometer. 



Results and Evaluations 

Table 1 shows the experimental conditions used and some results 
derived from them. 

Figs. 1 -3 show measured chlorine profiles and theoretical curves 
fitted to them. In type A experiments (PCVD deposited layers in a 
silica tube) the starting chlorine profile had to be taken into account 
when calculating the diffusion coefficient, because chlorine had al- 
ready diffused somewhat during the deposition time. Type B ex- 
periments (runs 4 and 7) allowed the solubility in equilibrium with 
the gas phase to be calculated, since it can be assumed that the 
innermost surface, which is always in contact with the chlorine 
atmosphere, should therefore be in equilibrium with it. 

Fig. 4 shows the temperature dependence of the diffusion coef- 
ficient of chlorine, in silica glass (HERALUX w.g. or deposited b> 
the PCVD process). The fitted straight line, showing no difference 
of the diffusion coefficient in the two different materials, corre- 
sponds to 

logjo (Dan" 1 s) = -11080 K/r- 3.498 . (-» 

The activation energy of Eq. (2) (212 kJ/mol) is rather similar io 
that estimated by Greeuw and Hasper [2] (193 kJ/mol), but our 
diffusion coefficients are about a factor 30 higher at the tempera- 
tures used by Greeuw and Hasper. Using the diffusion coefficient 
and the solubility of run 4 the total amount of chlorine expected 
to be taken up by the silica rod in the ampoule could be calculated. 
The expected value (23.4 ug) and that found (23 J ug) were in good 
agreement. 



Table 1 

Experimental conditions and some results derived from them. 
Type: (A) deposited PCVD layers in a tube or (B) elemental chlorine in silica ampoules or tubes. ?(C1 2 ): G 2 partial pressure during annealing. T: temperature, 
r: Annealing time at pressure P(Cl 2 ) and temperature T. S(C1 2 ): Equilibrium solubility of chlorine (weight per cent) in HERALUX w.g. at pressure PiO? 
and temperature T. S 0 (Cl 2 ): Equilibrium solubility of chlorine (weight per cent) at pressure 1 bar and temperature T. 0(C1 2 ): Diffusion coefficient. Analyst 
either SEM or neutron activation (NA). W (Cl 2 ): total weight of chlorine taken up by a silica rod of 2.40 cm 2 surface 
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Fig. 3 

Chlorine diffusion profile measured (symbols, run 7, see Table 1). 
The solid line has been fitted to the points using the diffusion 
coefficient calculated from Eq. (2) 



Chlorine diffusion profiles measured (symbols, runs 1-3, see 
Table 1) and calculated (solid lines fitted to the experimental 
points). 

O: T m 1965 K; x : T - 1748 K; A: T « 1628 K 
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The diffusion coefficient D of chlorine in silica glass as a function 
of the temperature 7" (runs 1 —4) 



Chlorine diffusion profiles measured (symbols, runs 4—6, Table 1). 

x : i>(Clj) = 13.4 bar, ■: P{C\ 2 ) - 5.6 bar. A: P{C\ 2 ) = 2.0 
bar. The solid line at 13.4 bar (•. x ) is fitted to the experimental 
points and the lines at 5.6 bar (■) and 2.0 bar (A) calculated from 
the first curve using the pressure dependence 



Runs 8-10 allow the Cl 2 pressure dependence of the solubility 
to be calculated. This is approximately proportional to the power 
of 1. 4 (Fig. 5). This dependence is in good agreement with runs 
4-6 given in Fig. 2, where the calculated curves according to that 
CI; pressure dependence are shown as solid lines. Such a Cl 2 pres- 
sure dependence is rather unexpected. It means that each Cl 2 mol- 
ecule generates four independent centres in the glass matrix in 
which it is dissolved. Two centres might be the chlorine atoms 
themselves, but the nature of the two additional centres Z is not 
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Fig. 5 

Cl 2 pressure dependence of the chlorine solubility in silica glass 
derived from runs 8—10 
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clear at the moment. From the reaction between gas phase and 
glass matrix 

Cl 2 = 2CI + 2Z (3) 
a mass action constant 

K m [Z] J [C1] 2 // > (C1 2 ) (4) 

is derived (brackets mean mole fractions in the glass). With 

CZ] - [CI] (5) 

the pressure dependence of the chlorine solubility is 

S(C1 2 ) - KP{C\ 2 ) 1 '* . (6) 

Using this relation the solubilities could be calculated at 1 bar CI 2 
pressure (Table 1). 

Conclusions 

Chlorine diffusion profiles in PCVD deposited silica glass 
and HERALUX w.g. agree quite well (see Fig. 4). On the 
other hand, the concentration of chlorine dissolved in 
PCVD silica is much higher than expected from the small 



Cl 2 pressure during deposition (1 -3 mbar). This might be 
a consequence of the very fast plasma reaction. It is well 
known that plasma reactions may have yields different from 
those expected for thermodynamic equilibrium [4], 

The pressure dependence of the solubility of chlorine in 
silica glass is rather strange. It should be an interesting in- 
vestigation to study the nature of the centres formed in silica 
when chlorine is dissolved in it. 

The authors are indebted to Mrs. B. Traftord for PCVD depo- 
sitions, to Mr. K. HerfT for sample preparation and to Mr. H. 
Jaspers for neutron activation analyses. 
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Additivity Rule of the Partial Molal Volume, 2. The Role of the Different 
Packing Efficiencies around Methylene and Methyl Groups of Normal 

Alkanes 

Y. Yoshimura and M. Nakahara 

Department of Chemistry, Faculty of Science, Kyoto University, Kyoto 606, Japan 

Crystals / Liquids I Mola! Volume / Solutions I Thermodynamics 

The molecular basis of the additivity rule of the partial molal volume at constant temperature and pressure has been studied for 
n-alkanes and it is found that the different packing efficiencies around the methylene and methyl groups play an important role in 
determining the partial molal volume. The different packing efficiencies have been examined in detail in the crystalline phases and the 
liquid phase. The methyl group is packed looser than the methylene group except for the triclinic crystal phase, and the temperature and 
pressure dependence of the packing efficiency of the methyl group is Jarger than that of the methylene group in the liquid state. The 
different packing efficiencies have been successfully interpreted in a semimicroscopic manner by taking account of the interfacial tension 

around the hard core of the solute molecule. 



L Introduction 

It is impossible to attain the firm and reliable compre- 
hension of the volume profile of the reaction in solution 
without knowing how the partial molal volume responses 
to a change in volume of the molecule dissolved in solution. 
The partial molal volume in solution or the density is one 
of the fundamental thermodynamic quantities to be deter- 
mined experimentally, while great efforts are being made to 
develop liquid theory so as to accomplish the ideal of de- 
riving the entire equation of state for any kind of molecules. 
Additivity is the simplest form as such a response and very 
easy to access empirically, and expected to contribute to the 
understanding of the volume profile of a homologous series 
of reactions in solution. The first paper [1] in this series has 
clarified the validity and limitations of the additivity rule 
of the partial molal volume for n-alkanes over a wide range 
of temperature and pressure. In the present paper we at- 



tempt to elucidate the physical significance of the parame- 
ters involved in the additivity rule from the molecular point 
of view by resorting to some model and crystal systems. 

In the previous work on n-alkanes [1] the following facts 
have been revealed: (i) in both solutions at infinite dilution 
and neat liquids, the additivity rule of the configurational 
partial molal volume at infinite dilution ( F? onf ) or the molal 
volume holds accurately in the form: 

v°-k t rt in 

= aw + P , l-> 

where F°, k t , rt, T, and /i denote the partial molal volume 
at infinite dilution (or the molal volume), the isothermal 
compressibility, the gas constant, the temperature, and the 
carbon number, respectively, and a and /? are the parameters 
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